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ABSTRACT 

In maximum-likelihood analyses of the Local Group (LG) acceleration, the object 
describing nonlinear effects is the coherence function (CF), i.e. the cross-correlation 
coefficient of the Fourier modes of the velocity and gravity fields. We study the CF 
both analytically, using perturbation theory, and numerically, using a hydro dynamic 
code. The dependence of the function on f2 m and the shape of the power spectrum 
is very weak. The only cosmological parameter that the CF is strongly sensitive to 
is the normalization as of the underlying density field. Perturbative approximation 
for the function turns out to be accurate as long as as is smaller than about 0.3. For 
higher normalizations we provide an analytical fit for the CF as a function of a% and 
the wavevector. The characteristic decoherence scale which our formula predicts is an 
order of magnitude smaller than that determined by Strauss et al. This implies that 
present likelihood constraints on cosmological parameters from analyses of the LG 
acceleration are significantly tighter than hitherto reported. 

Key words: methods: numerical, methods: analytical, cosmology: theory, dark mat- 
ter, large-scale structure of Universe 



1 INTRODUCTION 

The dipole anisotropy of the Cosmic Microwave Background 
(CMB) temperature is widely believed to reflect, via the 
Doppler shift, the motion of the Local Group (LG) with 
respect to the CMB rest frame. When transformed to the 
barycentre of the LG, this motion is towards (I, b) = (276° ± 
3°, 30° ±2°), and of amplitude «lg = 627 ±22 km-s" 1 , as in- 
ferred from the 4-year COBE data (Lineweaver et al. 1996) . 
Alternative models which assume that the dipole is due to 
a metric fluctuation (e.g., Paczyhski & Piran 1990) have 
problems with explaining its observed achromaticity and the 
relative smallness of the CMB quadrupole. 

An additional argument in favour of the kinematic in- 
terpretation of the CMB dipole is its remarkable alignment 
with the LG gravitational acceleration (or gravity) , inferred 
from galaxy distribution. The acceleration on the LG, in- 
ferred from the IRAS PSCz survey points only ~ 13° away 
from the CMB dipole apex (Schmoldt et al. 1999; hereafter 
S99, Rowan- Robinson et al. 2000). The alignment between 
the two vectors is expected in the linear regime of gravita- 
tional instability (Peebles 1980), and under the hypothesis 
of linear biasing between galaxies and mass. The ratio of 
the amplitudes of the velocity and gravity vectors is then 
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a measure of the quantity /3 — S7^ 6 /b, where Q m and b 
are the cosmological density of nonrelativistic matter and 
linear bias parameters, respectively. Therefore, comparisons 
between the LG gravity and the CMB dipole can serve not 
only as a test for the kinematic origin of the latter but also 
as a measure of j3. Combined with other constraints on bias, 
they may yield an estimate of Q m itself. 

However, the linear estimate of the LG velocity from a 
particular redshift survey will in general differ from its true 
velocity. The reasons are the finite volume of the survey, 
shot noise due to discrete sampling of the galaxy density 
field, redshift-space distortions and nonlinear effects. In a 
proper process of the LG velocity-gravity comparison, all 
these effects should be accounted for. 

A commonly applied method of constraining cosmolog- 
ical parameters by the LG velocity-gravity comparison is a 
maximum-likelihood analysis, elaborated by several authors 
(especially by Strauss et al. 1992, hereafter S92; see also 
Juszkiewicz, Vittorio & Wyse 1990, Lahav, Kaiser & Hoff- 
man 1990, S99). In this approach one maximizes the likeli- 
hood of particular values of cosmological parameters given 
the observed values of the LG velocity and gravitational ac- 
celeration. This enables one to constrain f3 and the relative 
amount of power on large scales within the framework of a 
given cosmology. 

The analysis of S92 constrained (3 to lie between 0.4 
and 0.85 (1 a). The acceleration on the LG was derived 
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there from the 1.2 Jy survey of IRAS galaxies. S99 repeated 
this analysis, with the LG gravity inferred from the recently 
completed IRAS PSCz catalogue. This catalog contains al- 
most three times more galaxies than its 1.2 Jy subsample. 
Still, the errorbars on 0, obtained by S99 (J3 = 0.70t°;| 5 at 1 
a), are not smaller than those obtained by S92. The volume 
surveyed is larger, shot noise is suppressed, but the errors 
remain big. Why? The authors blame nonlinear effects. 

In nonlinear regime non-local nature of gravity is un- 
veiled and the local relationship between the acceleration 
and velocity vectors is partly spoiled. In other words, the 
nonlinear velocity-gravity relation at a given point has scat- 
ter. As a result, the precision of determining f3 by the 
method described above is fundamentally limited, regard- 
less how well we can measure the LG gravity (and velocity). 

This argument sounds reasonable. However, if nonlinear 
effects are so strong, why is the misalignment angle between 
the LG gravity and velocity so small? Doesn't it actually 
suggest otherwise? This motivated us to reanalyze nonlinear 
effects in the LG velocity-gravity comparison. 

In a maximum-likelihood analysis, a proper object de- 
scribing nonlinear effects is the coherence function (here- 
after CF), i.e. the cross-correlation coefficient of the Fourier 
modes of the gravity and velocity fields (see S92 for details) ||] 
S92 devised a formula for the CF, calibrating it so as to fit 
the results of N-body simulations of a standard CDM cos- 
mology. They adopted this form of the function in all sub- 
sequent analyses, thus ignoring any possible dependence it 
may have on the shape of the power spectrum and its nor- 
malization. A similar approach was adopted by S99, who 
followed 'the S92 assumption that the CF does not change 
appreciably with the background cosmology'. S99 applied 
the same form of the function in two different cosmological 
models: spatially flat CDM cosmologies with respectively 
zero and non-zero (Q,a = 0.7) cosmological constant. The 
adopted value for the spectral parameter V was 0.25 and the 
spectra were cluster-normalized, so they had different values 
of as (r.m.s. mass fluctuations on the 8ft" 1 Mpc scale). 

By definition, on large enough, linear scales the CF is 
unity. On smaller scales we expect it to begin to depart 
from this value. The scale of departure marks a character- 
istic scale at which the non-local nature of gravity can no 
longer be ignored. Non-locality of the fields must be some- 
how correlated with their non-linearity, because it does not 
appear for linear fields. Since in more evolved models the 
nonlinearity scale is larger, it is natural to expect the CF to 
depend on the normalization of the underlying power spec- 
trum. If this is indeed the case, then this dependence should 
be modelled, to be accounted for in future LG velocity- 
gravity comparisons. The dependence on other cosmologi- 
cal parameters should also be studied. This is the aim of 
the present paper. It is organized as follows. In Section || 
we calculate the CF perturbatively. In Section |^ we describe 
the numerical simulations which we use to estimate the CF 
numerically. Results of the simulations are presented in Sec- 
tion ^. In particular, we compare the analytical estimates 
with the numerical estimates of the CF. In Section H we 
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show that the CF significantly affects the estimation of cos- 
mological parameters in maximum-likelihood analyses of the 
LG acceleration. Moreover, we constrain the CF in an alter- 
native way, adopted by S92. Summary and conclusions are 
in Section ^| 



2 ANALYTICAL CALCULATIONS 



The CF is defined as 

<gk • K) 



C(k) 



(i) 



(|gk| 2 )V 2 (|v k l 2 )V2 ' 

where gk and Vk are the Fourier components of the grav- 
ity and velocity fields, and (. . .) means the ensemble av- 
eraging. (Note that the definition of S92 lacks the com- 
plex conjugate sign.) The function can be interpreted as the 
cross-correlation coefficient of the Fourier modes of the grav- 
ity and velocity fields. It is important for the LG gravity- 
velocity comparisons, because it appears in the likelihood 
function for the LG velocity and acceleration (see Section ^J). 
As argued in the Introduction, the CF is also interesting in 
its own right - it carries information about non-locality, and 
indirectly about non-linearity, of the fields. In this section we 
calculate the CF for the fields which are weakly non-linear. 

By definition, C(k) is invariant to scaling of gk and Vk 
by an arbitrary constant. We are then free to choose the 
fields scaled so as to fulfill the following equations: 



Vg 
Vv 



5, 



(2) 
(3) 



Here, 8 is the mass density contrast and 9 is the velocity 
divergence, scaled in such a way that in the linear regime 
9 = 5. 

The gravity field is strictly irrotational. Hence 



(4) 



where <5k is the Fourier transform of the density contrast. 
Due to Kelvin's circulation theorem, the cosmic velocity field 
is vorticity-free as long as there is no shell crossing. Since ap- 
preciable shell crossing does not occur for weakly nonlinear 
fields, we have 



Vk = F ' 
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where #k is the Fourier transform of the velocity divergence. 
Using equations (^) and (^) we obtain 



C(k) 



(MO 
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<5k<5k} 1/2 <0k0k*} 1/2 ' 
We assume here that the density and velocity diver- 
gence fields are homogeneous and isotropic random fields. 
For such fields, the CF depends only on the magnitude of 
the wavevector. Furthermore, the density and velocity diver- 
gence are real fields, what implies that 5k* = <5-k, and sim- 
ilarly for the velocity divergence. All this implies that the 
CF is a real function, as follows. Namely, we have C(k) = 
C(-k), hence (* k 6» k *> = (5_ k 0_ k *) = (<5 k *0 k ) = (Mk*)*. 
Therefore, C(k) = [C(k)]*, that is C{k) is real. We may 
then cast it to an explicitly real form: 



C(k) = 



(fcgk) + (gkgk) 
2<Sk5k} 1 ' /2 <0k0k*} 1 / 2 



(7) 
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This formula is a good approximation to the CF as long 
as the velocity field has negligible vorticity. N-body simu- 
lations (Bertschinger & Dekel 1989, Mancinelli et al. 1994, 
Pichon & Bernardeau 1999) show that even in the case of 
fully nonlinear fields, the generated vorticity is small. 
We now expand 5 k and (9 k in perturbative series: Sk 



and similarly k = OP + # k ' + 8^' + 



,(2). 



g(3) 



8, we have (9 k = 5 k 



(i) 



.... Since in the linear regime 
We assume here that the initial (linear) density fluctuation 
field is a Gaussian random field. For such a field, all odd- 
order moments of the density contrast and of the velocity 
divergence vanish. Then, (<5 k 0k) = (si 1) 8^ ) *') + (8^ ) 8^ ) *) + 

(<W*) + (W) +0(a 6 ), where a 2 =(8 2 ), and 
similarly for other terms appearing in formula (m). Up to 
the leading order corrective term, this expansion yields 



C(k) = 1 - 



n(2)| 2 
k 



2(27T) a P(fc) 



(8) 



Here, P(fc) is the linear power spectrum, defined as 

(ar)-"<W>. 

The above formula is somewhat similar to that describ- 
ing weakly nonlinear corrections to the evolution of the 
power spectrum (Makino et al. 1992, Jain & Bertschinger 
1994). There is, however, also an important difference. 
Namely, all terms of the sort (a^ [3^ ), where a and /3 
stand for either 8 or 6, have remarkably cancelled out. In 
other words, unlike the weakly nonlinear power spectrum, 
the weakly nonlinear CF is constructed solely from second- 
order terms. 

To proceed further, we need the forms of <5 k 2 ' and (9 k 2 ' . 
Second-order solutions for the density contrast and (scaled) 
velocity divergence have been shown to depend extremely 
weakly on tt m and SIa (Bouchet et al. 1992, Bouchet et 
al. 1995, Bernardeau et al. 1995). This is also true for higher 
orders (see App. B.3 of Scoccimarro et al. 1998). Here we 
neglect the weak ^-dependence. Then (Goroff et al. 1986) 



and 

ei 2) - 



d 3 fcid 3 fc 2 
(2tt)3 



fo(ki + k 2 - k) J (2) (ki, ka) Sgsg 



(9) 



d 3 fcid 3 fe 



Mki + k 2 - k)K (2) (ki, ka) S^'S^ , (10) 



(l)r(l) 



(27T) 3 

where 8d is the Dirac delta 
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Hence we have 
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L(ki,k 2 ) 
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To determine (|<5 k 2) - 6> k 2) | 2 ), we need to evaluate the 

four-point correlations of the linear density field 8^ . For a 
Gaussian random field, 
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Using this property, it is a standard perturbative calculation 
to show that 



r(2) 
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d 3 gP(g)P(|k-q|)L 2 (q,k-q),(16) 



with L given by equation (h 

Note that j' 2 ' and K < - 2 >ha,ve first-order poles as ki — > 
or k 2 -> for fixed k: J (2) ~ AT (2) ~ (1/2) cos0(fci/fa + 
k 2 /ki), where $ is the angle between ki and k 2 . However, 
in the expression (|l^) for 5 k 2 ' — 8^ they cancel out: the func- 
tion L has no poles. This results in significant simplification 
of the integrand in equation (jl6|) . 

Using formula (H) we obtain 



C(k) 



-p-\k) 



d 3 q 
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-P(g)P(|k-q|)L 2 (q,k-q).(17) 



Now we write the integral in spherical coordinates q, i? and 
4>, the magnitude, polar angle, and azimuthal angle, respec- 
tively, of the wavevector q. Then with the external wavevec- 
tor k aligned along the z-axis the integral over (f> is trivial 
and simplifies J d 3 q to the form 2-7T J dqq 2 J d cos '8. Further- 
more, we have 



L 2 (q,k-q) 



49 [l + (g/fc)2-2( 9 /fc) M ] : 



(18) 



where fj, — cos i9. This suggests a change of variables s = q/k. 
Performing this finally yields 



1 - C{k) 



with 
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49(2^) 2 P(fc) Jkmin/h 
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/k pi 
dss 2 P(ks) / d/iX 



(1 + s 2 -2s(j,y 



(19) 



(20) 



The quantities k m i n and k max are the cutoffs of the power 
spectrum. Physically, they are related to the effective depth 
of the survey, from which the specific spectrum is extracted, 
and the virialization scale, respectively. (For k > k max per- 
turbative expansion breaks down.) In numerical simulations, 
kmin and k m ax are determined by the simulation box size 
and elemental cell size, respectively. Details will be given in 
the next Section. 

Given the linear spectrum P(k), equation (jl^) provides 
the second-order correction to the linear value of the CF. 
Results for the PSCz spectrum and for the standard-CDM 
spectrum are presented in Section ^| 



4 Chodorowski & Ciecielg,g 



Table 1. Parameters of the runs displayed in Figure [lj. The 
acronym l-l stands for 'low N A , low fejsrq', h-l stands for 'high 
N 3 , low /cNq'i an d h-h stands for 'high N A , high fcjjq'- 



run 


grid 


box size [h 1 Mpc] 


fc Nq [hMpc- 1 ] 


l-l 


64 3 


50 


4. 


h-l 


128 3 


100 


4. 


h-h 


128 3 


50 


8. 



3 NUMERICAL SIMULATIONS 

Following Peebles (1987), instead of using a N-body scheme 
we model cold dark matter as a pressureless cosmic fluid. 
Using the Eulerian code CPPA (Cosmological Pressureless 
Parabolic Advection, see Kudlicki, Plewa & Rozyczka 1996, 
Kudlicki et al. 2000, Kudlicki et al. 2001b) we solve its dy- 
namical equations on a uniform grid fixed in comoving co- 
ordinates. The main improvements of CPPA over the origi- 
nal Peebles' code are parabolic density and velocity profiles, 
variable timestep, periodic boundary conditions and a flux 
interchange procedure, implemented as an approximation to 
the solution of the Boltzmann equation. 

We chose to use a grid-based code rather than a N- 
body code because it directly produces a volume-weighted 
velocity field. This is important because in the definition of 
the CF, equation ([!]), the velocity field is volume- weighted, 
not mass-weighted. Moreover, the field is evenly sampled, 
which is convenient for FFT techniques. 

The linear velocity depends on the cosmological con- 
stant (Qa) very weakly (e.g. Lahav et al. 1991); this also 
holds for higher orders (see Bouchet et al. 1995, Appendix 
B.3 of Scoccimarro et al. 1998 and Nusser & Colberg 1998). 
Therefore, it was a good approximation to assume SIa = 
in our models. We have thus studied two zero-A models with 
Q m = 1 and Q, n = 0.3, assuming Gaussian initial conditions. 
The parameters of the runs are given in Table |l|. 

To make the simulated gravitational field as close as 
possible to that inferred from the IRAS PSCz survey, the 
mass power spectrum that we adopted was that estimated 
for the PSCz galaxies (Sutherland et al. 1999): 

Bk 



p(k) = 



(21) 



{1 + [ak + {bkf/ 2 + (cfc) 2 ]"} 2 /" 
a = 6.4/r/!, _1 Mpc, fe = 3.0/r/!, _1 Mpc, 
c = 1.7/r Ji _1 Mpc , ^ = 1.13, 

with r = 0.2 as best fitted value. The power spectrum em- 
ployed in simulations is effectively truncated at both large 
(corresponding to the lower cutoff k m in) and small (cor- 
responding to the upper cutoff k max ) scales. Specifically, 
k m in = 2-k/L, where L is the simulation box size, and 
k m ax = &N q = (N/2)k m in- Here, A?Nq is the so-called Nyquist 
wavevector and JV 3 is the grid size. 

To normalize the power spectrum we used the observed 
local abundance of galaxy clusters. The present value of erg, 
labelled os,o, i s a function of O m and for the case of Qa = 
it is estimated by the relation (Eke, Cole & Frenk 1996) 

(0.52 ± 0.04)n-°- 46+0 - lon '" . 




O"8,0 



(22) 



This relation changes only slightly with the shape of the 



[h/Mpc] 



Figure 1. The CF for three values of as, for three Einstein— de 
Sitter simulations with different resolutions. (Their parameters 
are given in Table |lj) Dotted lines show the results of run l-l, 
dashed lines h-l, and solid lines h-h. The CF is similar for runs 
l-l and h-l. For k < 2 ftMpc -1 , it is similar for all runs. For each 
run, the function bends down at k ~ 0.5&Nq- 



power spectrum. It is also very similar for the case of non- 
zero S1a, flat models (Ha = 1 — fim)- For = 1, as,o — 
0.52, while for Q m = 0.3, cr 8 ,o - 0.87. 



4 RESULTS 

Firstly, we tested the dependence of the results on resolu- 
tion. For this purpose we performed three simulations of an 
Q m = 1 Universe with the PSCz power spectrum (see Ta- 
ble ^). The runs l-l and h-l have the same spatial resolution 
but different grid resolutions, while the run h-h has higher 
spatial resolution. In Figure we present the temporal evo- 
lution of the CF. Specifically, we show it for 3 values of as: 
0.1, 0.3, and 0.52. As expected, the results of runs l-l and h-l 
do not differ significantly. In contrast to the run h-l, the run 
l-l has no modes corresponding to the scales greater than 
50 h^ 1 Mpc. These scales, however, are well in the linear 
regime, so C = 1 to good accuracy. 

In grid simulations, the largest wavevector is the 
Nyquist wavevector, corresponding to the Nyquist wave- 
length, i.e., the smallest wavelength, of two cells. Its values 
for the three runs are given in the last column of Table JjJ 
Fourier modes with k > fcNq do not have physical meaning. 
Inspection of Figure [j], however, shows that the CF bends 
down already at k ~ 0.5/cNq- Such a scaling of the bending 
point with the Nyquist wavevector strongly suggests that 
this is a resolution effect. Given the similarity of the three 
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k [h/Mpc] 

Figure 2. The CF for simulations 128 3 with the box-size of 50 
h Mpc. Dotted lines are for Q m = 1.0, dashed for f2 m = 0.3, 
and solid show predictions of perturbation theory for Q m = 1.0, 
according to formula dlj|). The CF depends on Q m very weakly. 
For <xg < 0.3, it is well predicted by the second-order approxima- 
tion, except for the highest values of the wavevector. 



curves for k < 2 /iMpc -1 , we can expect the CF to be free 
of resolution effects up to k ~ 0.5&N q , for any grid. 

Next, we test the perturbative approximation for the 
CF. In Figure |^, we show the CF from simulations with the 
grid 128 3 and the box-size of 50 /i -1 Mpc, for wavevectors 
up to 0.5fcN q = 4.0 /iMpc" 1 . Dotted lines are for f2 m = 1.0, 
dashed for fl m — 0.3, while solid show predictions of per- 
turbation theory for Q. m = 1.0, according to formula (|l9|). 
(We numerically evaluate the integral in this formula.) The 
CF depends on fl m very weakly. Moreover, except for the 
highest k- values, it is well predicted by the second-order ap- 
proximation, as long as erg is smaller than ~ 0.3. We have 
checked that for erg < 0.3 the quantity 1 — C(k) scales ap- 
proximately like as, as predicted by equation ( |l9| ) (the nor- 
malization of P(k) is proportional to erf). For higher values 
of as, the second-order approximation overestimates the de- 
viation of C from the linear value, unity. This behaviour is 
similar to the nonlinear evolution of the power spectrum, 
overestimated by second-order terms (Jain & Bertschinger 
1994). 

To describe the CF for erg > 0.3, we analyzed it for 
23 output times of the high-resolution simulation h-h, cor- 
responding to the values of erg in the range 0.3 < as < 1.0. 
We found that the erg-dependence of C can be well fitted as 

C(fc) = exp (-ak) , (23) 

where 




k [h/Mpc] 



Figure 3. The coherence function for cluster-normalized cos- 
mological models. Dotted lines show the function from the nu- 
merical simulation h-h, while solid ones are the fits according to 
formula 



f 0.757 af for as < 0.3, 

\ -0.059 + 0.423 a s for 0.3 < a s < 1.0 . 1 ' 

In Figure |^ we show the CF from the simulation h-h, for erg 
given by the cluster normalization. (Having shown very weak 
dependence of the function on f2 m , in this plot we present 
the results of Einstein-de Sitter simulations only, which are 
simpler to evolve numerically.) In particular, in this figure 
the dotted line for erg = 0.52 corresponds to the lowest solid 
line in Figure The fits are good; we have checked that for 
other values of erg they are as good as those shown here. 

To test the dependence of the CF on the underlying 
power spectrum, we performed a 64 3 simulation with the 
standard CDM spectrum (r = 0.5). The box-size was 50 

Mpc, so 0.5 fcNq = 2. The results are shown in Figure ^j. 
Solid lines are drawn according to our fit (eq. ||), obtained 
for the PSCz spectrum, while dotted lines are from the sim- 
ulation. A good agreement between them apparent in the 
figure implies that the dependence of the CF on the power 
spectrum is very weak, at least for a CDM-like family of the 
spectra. 

The CF has been modelled by S92, who calibrated it 
so as to fit the results of N-body simulations of a standard 
CDM cosmology. In Figure |E| we show S92's prediction for 
C, as well as our predictions, for the standard CDM power 
spectrum and erg normalization of S92 (0.625). The discrep- 
ancy of our results with the formula of S92 is drastic! Instead 
of a characteristic decoherence scale of 4.5 ft" 1 Mpc (S92), 
our formula ( pj|) suggests a fraction of a megaparsec. We 
will comment on that in the summary. 
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Figure 4. The coherence function from a simulation with the 
standard-CDM power spectrum (dotted) versus our formula (Bq), 
obtained for the PSCz spectrum (solid). 
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Figure 5. The coherence function for a standard-CDM cosmo- 
logical model with erg = 0.625. The dotted line shows the function 
from our standard-CDM simulation, while the solid one is the fit 
according to formula (Kit). The dashed line is the formula (18) of 
S92, with r c = 4.5 h'^Mpc. 



5 FROM COHERENCE TO PROBABILITY 
CONTOURS 

The aim of this Section is twofold. Firstly, we show that the 
CF is very important in the analyses of the LG acceleration, 
because, together with other factors, it determines the rel- 
ative likelihood of different cosmological models. Secondly, 
we show that the probability distribution for the LG ac- 
celeration amplitude and the misalignment angle, resulting 
from our estimate of the CF, is consistent with that obtained 
from mock IRAS catalogs, constructed by S92. In contrast, 
adopting S92's formula for the function results in a distribu- 
tion which is too broad. Here we only outline the necessary 
formalism; for more details the reader is referred to S92 and 
Chodorowski & Ciecielag (2001). 

Let /(g, v) denote the probability density distribution 
for observing particular values of the LG gravity and ve- 
locity, given some assumed CDM cosmological model. The 
model is fully specified by the value of /3, the power spec- 
trum shape parameter F, and the normalization as. Due to 
Bayes theorem, it is possible to evaluate from / the relative 
likelihood, C, of different models (/3, T, as). The likelihood 
function is usually defined as 



£ = /■ 



(25) 



Here, the observed values of the Local Group velocity and 
gravity are inserted in / and serve as constraints on cosmo- 
logical parameters. 

As a functional form of /, S92 and S99 adopt a multi- 
variate Gaussian. This assumption has support from numer- 
ical simulations (Koffman et al. 1994, Kudlicki et al. 2001a), 
where the measured nongaussianity of g and v is small. This 
is rather natural to expect since, e.g., gravity is an integral 
of density over effectively a large volume, so the central limit 
theorem can at least partly be applicable (but see Catelan 
& Moscardini 1994). However, the approximate Gaussianity 
of g and v by no means implies that the fields are linear. In 
contrast, the above considerations suggest a Gaussian ap- 
proximation for the form of /, but with the covariance ma- 
trix calculated accounting for the nonlinear effects. 

After some algebra (see e.g. Juszkiewicz et al. 1990), 
the joint distribution for g and v can be cast to the form: 



/(&v) = 



(1 



r 2 ) 



-3/2 



(27T) 



3n-?n-3 



■ exp 



x + y — 2rfixy 



2(1 



(26) 



where a g and a 



denote the r.m.s. values of a single spatial 
component of gravity and velocity, respectively. From sta- 
tistical isotropy, a 2 g — (g ■ g)/3, and a% = (v ■ v)/3. Next, 
(x,y) = (g/cr 9 , v/<Tu), and fi = cos (9 with 9 being the mis- 
alignment angle between g and v. Finally, r is the cross- 
correlation coefficient of gi with Vi, where g% (vi) denotes an 
arbitrary spatial component of g (v). From isotropy, 



(g' v ) 



(g 2 ) 1 /2(„2)l/2 • 

Also from isotropy, 

(Xiyj) = r5 i:j , 



(27) 



(28) 



where Sij denotes the Kronecker delta. In other words, there 
are no cross-correlations between different spatial compo- 
nents. 
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In the limit of linear fluctuations and with perfect sam- 
pling of the density field, the distribution (|2^) reduces to 
<5d(x — y)A/3(v), where A/3 denotes a tri-variate normal dis- 
tribution. In a real world, there are nonlinear effects (NL), 
as well as finite- volume effects (FV), which make the cross- 
correlation coefficient deviate from unity. In a related paper 
(Chodorowski & Ciecielag 2001) we will show that r can be 
approximately written as 



r = r N L ■ r F v • 



(29) 



i.e. that the contributions to r from nonlinear effects sepa- 
rate from those from finite volume. In the present paper we 
study nonlinear effects. The cross-correlation coefficient due 
to them is 



r o °° c(k)w g (k)w v (k)p(k) dk 

/ °° W g {k)W v (k)P{k) dk 



(30) 



Here W g and W v are the observational filters of g and v 
(for details see S92). Thus, nonlinear effects enter into the 
correlation coefficient of the joint distribution function via 
the CF. In the linear regime C = 1 and, if we neglect other 
effects, r = 1. Since the distribution function is directly re- 
lated to the likelihood of different world models (eq. Em, in 
analyses of the LG acceleration, the CF affects the estima- 
tion of cosmological parameters. 

In equation ( |3o| ) the CF is multiplied by the windows 
through which the gravity and velocity of the LG are mea- 
sured. Therefore, smoothing effectively filters out the high-fc 
tail of the coherence function. It is instructive to write 



/ °° [I - C(k)] W g (k)W v (k)P{k) dk 
$™W g {k)Wv{k)P{k)&k 



(31) 



the stronger the deviation of C from unity so is the de- 
viation of tnl- In Figure ^| we plot the integrand [I — 
C(k)]W g (k)W v (k)P(k) as a function of k. (Strictly speak- 
ing, we plot the function [I - C(k)]W g {k)W v (k)kP(k), be- 
cause the fc-axis is logarithmic.) Here C(k) is given by our 
fit (eq. ^), P(k) is the spectrum of the PSCz galaxies, and 
the velocity window is that introduced by S92, with a small- 
scale cutoff, r m in = 1 ft -1 Mpc, to reflect the finite size of the 
LG. The gravity window is either the standard IRAS win- 
dow with a small-scale smoothing r s = 5 h~ l Mpc (S92), 
or equal to the velocity window. The first case corresponds 
to the present situation, where the inferred LG gravity is 
commonly smoothed with the standard IRAS window, and 
results in the dashed curve. The second case describes an 
ideal situation, where the galaxy distribution around the 
LG is sampled so densely that there is no need to smooth 
its gravity beyond the size of the LG. That case results in 
the dotted curve. We see from the figure that at present 
it is sufficient to know the behaviour of the CF up to at 
most k = 2 ft Mpc -1 , and it will never become necessary 
to know it for k > 4 /iMpc -1 . This is why we have set 
up the resolution of the h-h simulation in such a way that 
0.5fcNq = 4 /iMpc" 1 (Section §). 

In contrast to our approach, S92 did not determine 
the CF from its definition. Instead, using standard CDM 
N-body simulations, they created mock IRAS catalogs for 
'observers' (N-body points) selected with similar properties 
to those of the LG (primarily the velocity). For each of these 
observers, S92 computed the IRAS acceleration of the LG, 
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Figure 6. The integrand in the numerator of the right-hand side 
of equation ( pl| ) as a function of the wavevector. Since the fc-axis is 
logarithmic, the ordinate is multiplied by and extra power of k, so 
equal areas under the function correspond to equal contributions 
to the deviation of the correlation coefficient from unity. Units of 
the ordinate axis are arbitrary. 



as the filter W g using the so-called standard IRAS window. 
They then matched the resulting distribution for the am- 
plitude of g and the misalignment angle with the probabil- 
ity contours resulting from a formula following from equa- 
tion (^fi|). Since the LG velocity was considered by S92 as a 
constraint, we need the conditional probability density func- 
tion /(g|v). It readily results from formula (|2q): 



/(g|v) = (2tt) 



-3/2 



'(I 



-3/2 



exp 



(x - ry) 2 
2(1 -r 2 ) 



(32) 



(Juszkiewicz et al. 1990). The distribution for g and 9 results 
from the above formula by multiplying it by 2irg 2 sin 9. 

The approach adopted by S92 implies that the finite 
volume effects are also present, so they should also be mod- 
elled. We did so, and the results are presented in Figures ^ 
and ||. Specifically, in both figures the scattered points are 
the distribution for g and 9 for the LG observers, simulated 
by S92 (Strauss, private communication). On this distribu- 
tion, we superimpose the probability contours resulting from 
equation (^), for observers constrained by the CMB dipole 
(v = komb). The contours are drawn with the finite- volume 
effects modelled according to S92. Figure |^ shows the con- 
tours drawn with the CF given by our formula (|2^), while 
Figure ^ shows them for C given by the fit of S92. We see 
that the fit of S92 results in obviously too broad probabil- 
ity contours. In contrast, our formula results in the contours 
which at first look seem to be consistent with the simu- 
lated distribution. At a closer inspection, one might worry 
that they are also (slightly) too broad: e.g., outside the 95% 
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Figure 7. Contours of probability density as a function of am- 
plitude of the LG gravity and the misalignment angle. The LG 
velocity (CMB dipole) is used here as a constraint. The contours 
correspond to the probability levels of 68%, 90%. and 95%. The 
coherence function is given by our formula (B3|). The scattered 
points in this Figure and in Figure |s| are S92's values from simu- 
lations. 



tightens up the contours slightly' (S92). Still, as very little, 
if any,0 modification is needed, it may be just enough. 

The above uncertainty has little relevance, since the 
aim of this section was qualitative rather than quantitative. 
We showed that the CF is important in analyses of the LG 
acceleration, and that an alternative way of calibrating it, 
adopted by S92, also points towards much smaller decoher- 
ence scale. We did not, however, attempt to determine the 
CF in this way. Actually, we think that such a determination 
is non-trivial, since one has to separate carefully the influ- 
ence on the correlation coefficient of nonlinear effects from 
remaining effects. 

In calculating Figures]^ and ^, a correction was adopted. 
We computed a 2 , and a v according to the linear theory 
[(67r 2 ) -1 J W 2 (k)P(k)dk]. Since on small scales the gravity 
window smoothes more heavily than the velocity window 
(S92), the resulting a g was smaller than a v . This had the 
effect that the predicted and simulated distributions were 
slightly off-set horizontally. To correct this, we equated a g 
to a v , obtaining Figures |^ and ^| 

Is this slight correction valid? Nonlinear gravity is 
known to be slightly larger than nonlinear velocity smoothed 
with the same filter (Berlind et al. 2000, Kudlicki et 
al. 2001a). Therefore it may well be that the effect of differ- 
ent filters compensates here, at least partly, with the nonlin- 
ear effect. This means that in the LG velocity-gravity com- 
parison there are remaining residual nonlinear effects, that 
still deserve further modelling. We plan to do this elsewhere. 




200 400 600 800 1000 
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Figure 8. As in Fig. ^, but with the coherence function given 
by the fit of S92. 



probability contour there are only 5 points out of 200. How- 
ever, the simulated distribution was constructed by S92 un- 
der an additional constraint of small shear of the velocity 
field around the LG, and our model does not include this. 
The effect of the local shear constraint 'is minor and only 



6 SUMMARY AND CONCLUSIONS 

We have studied the coherence function (CF), i.e., the cross- 
correlation coefficient of the Fourier modes of the cosmic ve- 
locity and gravity fields. This function describes nonlinear 
effects in maximum-likelihood analyses of the LG accelera- 
tion. It is important, since it affects the likelihood function 
and hence the estimation of cosmological parameters (Sec- 
tion ^). We have determined the CF both analytically using 
perturbation theory (Section^), and numerically using a hy- 
drodynamic code (Section ^|and ^). The dependence of the 
function on Q. m and on the shape of the power spectrum has 
turned out to be very weak. The only cosmological parame- 
ter that the CF is strongly sensitive to is the normalization, 
08, of the underlying density field. We have found that the 
perturbative approximation for the function is accurate as 
long as <T8 is smaller than about 0.3. For higher normal- 
izations we have provided an analytical fit for the CF as a 
function of ag and the wavevector. The characteristic de- 
coherence scale which our formula predicts is an order of 
magnitude smaller than that found by S92. 

To analyze the above discrepancy we have followed the 
approach of constraining the CF adopted by S92. Specifi- 
cally, we have calculated the distribution function for the 
amplitude of the acceleration of the LG and the misalign- 
ment angle, given the value of the LG velocity, and compared 
with that obtained from mock IRAS catalogs. The distribu- 
tion resulting from our estimate of the CF turned out to 



2 The discrepancy may be simply statistically insignificant. 
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be consistent with the simulated one. In contrast, adopt- 
ing S92's formula for the function resulted in a distribution 
which is too broad. We believe therefore that we have de- 
rived the correct form of the CF. The origin of the error in 
the analysis of S92 remains unclear to us. 

Tighter probability contours for the LG gravity imply 
tighter confidence intervals for estimated cosmological pa- 
rameters. The likelihood contours for the parameters can 
only be drawn given the data, thus we leave it for fu- 
ture analyses of the observed LG acceleration. This paper 
strongly suggests that with proper account for nonlinear ef- 
fects in such analyses, the value of j3 can be determined with 
significantly greater precisionF] than is currently believed. 
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